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SEMIEMPIRICAL ELECTRIC DIPOLE MOMENT AND ITS 
DERIVATIVE IN H 2 -H 2 AND H 2 -H COLLISIONS 
by R. W. Patch 
Lewis Research Center 

SUMMARY 

In a gas, when two Hg molecules or an Hg molecule and an H atom collide, a mo- 
mentary electric dipole moment results. This electric dipole moment may cause 
pressure -induced infrared absorption. A knowledge of the strength of this absorption is 
necessary for radiant-heat-transfer calculations at high pressures. The strength de- 
pends, among other things, on the value of the electric dipole moment and its partial 
derivative with respect to Hg internuclear distance. These quantities cannot be deter- 
mined accurately from ab initio calculations on present day computers if the intermolec- 
ular distance is much greater than 4 bohr (2. 117x1 0 - ^ m), nor can they be determined 
unambiguously from experiments. 

In this report the electric dipole moment and its derivative for Hg-Hg and H2-H 
collisions were calculated for intermolecular distances of 4 to 8 bohr (2. 117x1 O - *® to 
4. 233x10"^® m) by using a semiempirical model for Hg-Hg previously developed by 
others. In this model the dipole moment and its derivative are the sum of an overlap 
contribution and a quadrupole- induced contribution. The quadrupole- induced contribu- 
tion is analytically simple but requires correction of previously published expressions. 
The overlap contribution is analytically complicated. For Hg-Hg the overlap analysis 
was improved, and for H^-H it was reformulated. For both contributions the only ex- 
perimental input was the equilibrium internuclear distance of Hg, which was obtained 
from spectroscopic data. 

Dipole moment and its derivative were calculated for a number of collision configu- 
rations. For Hg-Hg four planar configurations and one nonplanar configuration were 
used. For Hg-H three configurations were employed. 

The semiempirical theory used in this report is not valid for intermolecular dis- 
tances less than 4 bohr (2. 117x10”^ m) because of neglect of configuration interaction. 
Therefore, to obtain a comparison with Hg-Hg experiment, results from this report for 
large intermolecular distances were faired into ab initio results for small intermolecular 
distances. The resulting integrated absorption coefficient at 298 K for the pressure- 
induced fundamental vibrational transition was 87 percent of the experimental value. 



INTRODUCTION 


In certain high-temperature propulsion devices including gas-core nuclear rockets, 
an important mechanism of heat transfer is radiant energy exchange. In the gas-core 
rocket this radiation occurs between the uranium plasma and the hydrogen gas and be- 
tween the hydrogen and the wall (refs. 1 and 2). Because of the high pressures in such 
devices, it is necessary to know the strength of pressure-induced infrared absorption 
of hydrogen to calculate the heat transfer (ref. 3). This strength is also necessary in 
constructing models of late-type stars (ref. 4). The strength depends on the interaction 
energies and the electric dipole moments occurring in the collisions of Hg molecules 
with various collision partners. In the cases mentioned, important collision partners 
include Hg and He and probably H. No experimental data taken much above room tem- 
perature exist for collision partners Hg and He, and none at all for H. The present in- 
vestigation is restricted to calculating the electric dipole moment and its derivative for 
collision partners Hg and H, and is part of a project to calculate pressure- induced vi- 
brational strengths for temperatures of 600 to 7000 K. 

Pressure-induced absorption by an Hg molecule occurs when a photon is absorbed 
because of the electric dipole moment resulting from the proximity of a collision part- 
ner. To conserve energy, the relative kinetic energy of the colliding molecules or the 
internal energy of Hg, or both, must change. If only the relative kinetic energy changes, 
the process is called pressure-induced translational absorption. If the Hg rotational 
energy changes but not the vibrational energy, the process is termed pressure- induced 
rotational absorption. If the Hg vibrational energy changes, the process is referred to 
as pressure-induced vibrational absorption. Vibrational absorption is the most impor- 
tant of these three processes in heat transfer at the temperatures in propulsion devices 
because it occurs at the shortest infrared wavelengths, which are nearest the peak of 
the Planck (blackbody) function. 

For large intermolecular distances the electric dipole moment which occurs when 
Hg collides with an atom or homonuclear diatomic molecule may be assumed to be the 
sum of two contributions. The quadrupole- induced contribution is due to the dipole mo- 
ment induced in the atom or in the homonuclear diatomic molecule by the quadrupole 
moment of Hg. If the collision partner is a homonuclear diatomic molecule, its quad- 
rupole moment may also induce a dipole moment in Hg. In this case, the two dipole 
moments are additive. The other contribution to the dipole moment comes from the 
distortion of the electron charge clouds of the two colliding molecules which results 
from their overlap. This contribution is much more difficult to calculate than the first 
and is the principal subject of this report. The overlap contribution is most important 
in pressure-induced vibrational absorption, which is the process emphasized herein. 

In vibrational absorption the derivative of the dipole moment with respect to Hg inter- 
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nuclear distance is of much greater importance than the value of the dipole moment it- 
self. 

A number of theoretical investigations have been relevant to pressure- induced vi- 
brational absorption by Hg. Van Kranendonk and Bird (ref. 5) gave a semiempirical 
theory for the overlap contribution in Hg-Hg collisions based on the distortion of the 
orbitals of two isolated repelling H atoms. They calculated what they called an 
"indirect” effect in the overlap contribution. When this was added to the quadrupole- 
induced contribution, their vibrational absorption coefficients were lower than measured. 
Britton and Crawford (ref. 6) added a "direct" effect to Van Krandendonk and Bird's 
results to get good agreement with absorption experiments at room temperature. 

Van Kranendonk (refs. 7 and 8) took quantum corrections to the pair distribution function 
into account to obtain an absorption coefficient theory applicable down to 1 50 K. None of 
these investigators considered configuration interaction because they were interested 
only in intermolecular distances of 4 bohr (2. 117x10”^ m) and greater. Kolos and 
Wolniewicz have since calculated more accurate values of Hg quadrupole moment and 
polarizability (refs. 9 and 10, respectively). Patch (ref. 11) has done an ab initio cal- 
culation of interaction energy and dipole moment and its derivative for Hg-Hg, which in- 
cludes configuration interaction and therefore is believed to be valid for intermolecular 
distances of 4 bohr (2. 117x10” ^ m) or less. These small distances become important 
well above room temperature. No investigations of pressure- induced Hg-H absorption 
have been made. 

This report calculates dipole moments and their derivatives for H 0 -H 0 and H„-H 

_ i o « « ^ 

for intermolecular distances of 4 bohr (2. 117x10 m) and greater. At such distances, 
accurate ab initio calculations are impractical on present day computers, so the semi- 
empirical model of references 5 and 6 was used. It was refined for Hg-Hg and had to 
be reformulated for Hg-H. 

The main body of the report is divided into two parts: one for Hg-Hg and for Hg-H. 
Each part is divided into an analysis section and a results and discussion section. The 
analysis section contains the derivation and points out the differences between it and 
references 5 and 6. The results and discussion section gives numerical results, com- 
parisons with previous Hg-Hg investigations, limitations, and potential applications. 

DIPOLE MOMENT AND ITS DERIVATIVE IN H 2 "H 2 COLLISIONS 


Analysis 

In the Born-Oppenheimer approximation and for intermolecular distances of 4 bohr 
(2. 117x10”^ m) or greater, the dipole moment of Hg-Hg and its derivative consist of 
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two additive contributions: an overlap contribution and a quadrupole- induced contribu- 
tion (refs. 5 and 6). In this section the overlap contribution is calculated from the 
semiempirical model of references 5 and 6 with the following refinements: 

(1) Previously neglected terms are included. 

(2) The repulsive distortion parameter X is based on a repulsive H-H calculation 

using the same orbital exponent as for Hg. 

(3) A nonplanar configuration is included. 

(4) The dipole moment, as well as its derivative, is obtained. 

The overlap contribution is then expanded in spherical harmonics. Corrected equa- 
tions for the expansion coefficients of the quadrupole- induced contribution to the deriva- 
tive of the dipole moment (including anisotropy of Hg polarizability) are given, together 
with expansion coefficients for the quadrupole- induced contribution to the dipole moment 
itself. All equations are in atomic units (bohr, hartree, and electron charge). 

Coordinates . - The coordinate system is shown in figure 1. The origin of the 
Cartesian coordinates x^, x^, and x^ is on the intermolecular axis halfway between 



Figure 1 . - Coordinates in H 2 -H 2 collisions. The molecules are a-b and c-d; protons are at a, b, 
c, and d ; the X3 axis passes through the midpoints of the molecules; 0^ and 0 2 are polar angles; 
(p^ and (p2 are azimuthal angles measured in the xj-x 2 plane. 


molecules a-b and c-d (symbols are given in the appendix). The relative positions of 
the molecules are given by polar angles angles 9^ and anc * azimuthal angles <p ^ 
and q> 2 ' Only the x^ component of the dipole moment and its derivative are considered 
for Hg-Hg because Britton and Crawford (ref. 6) showed that the other components have 
only minor importance. 

Overlap contribution to dipole m oment. - For intermolecular distances of 4 bohr 
(2. 117xl0~ 10 m) or greater, it is reasonable to neglect configuration interaction, so the 
antisymmetric system wave function 4/ was assumed to be 
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¥ = “ ^3 “ ^4 + ^5 


( 1 ) 


where the determinantal wave functions and i//g to are given by 

= (4I) _1/2 det[k(l)n(l),b(2)i3(2) ) c(3)«(3),d(4)^(4)]" 

= (4! ) -1 / 2 det[a(l)o'(l),b(2)^(2),c(3)i3(3),d(4)a(4)] 

► 

</^ 4 s (4!)" 1 /2 detfa(l)/3(l), b(2)a(2), c(3)a(3), d(4)/3(4)J 
= (4! ) -1 / 2 det|a(l)/3(l), b(2)ar(2), c(3)/3(3), d(4)a(4)]^ 


Here a(l) is an atomic orbital for electron 1 centered on proton a and similarly for 
b(2), etc. Also, «(1) is the spin eigenfunction of electron 1 with the component of spin 
angular momentum along the axis of quantization equal to h/2, and £(3) is the spin 
eigenfunction of electron 3 with the component of spin angular momentum along the axis 
of quantization equal to -h/2. 

The Xg component of the overlap dipole moment is given by 



where Xg^ is the Xg coordinate of proton p, and the p under the summation sign in- 
dicates summation over all protons. Similarly, Xg^ is the Xg coordinate of electron i. 
The integration is over configuration and spin space of the four electrons. Equation (3) 
readily simplifies to 


M 


ax 3 



x 3 i ^ dV 2 dV 3 dV 4 


J * 2 dVj dV 2 dVg dV 4 


( 4 ) 
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Equation (4) may be simplified by noting the properties of products of determinantal 
wave functions (ref. 12) and by introducing the following five types of integrals: 


®ab s j a (l) b OO dv i ( a and b in s ame molecule) 

(5) 

s ac s J' a(l)c(l)dvj (a and c in different molecules) 

(6) 

m 3aa s /' a2 (l) x 3al dv l 

(7) 

m 3ab S J' ^ l)b(l) x 3abl dv l ( a and 13 in same molecule) 

(8) 

m 3ac S J' a (l) c (l) x 3 ac i dv l ( a and c * n different molecules) 

(9) 

where x g a j is the Xg component of a vector from proton a to electron 1, x g ac j is the Xg 
component of a vector from the center of line ac to electron 1, and dvj is an element 
of volume for electron 1. In addition, it was assumed that 

f a 2 (l)d Vl = 1 

(10) 


and likewise for similar integrals. Hence, equation (4) becomes 



N 1 + N 2 + N 3 
E 


( 11 ) 


where 


N i s 


-4 (m 


3aa 


+ m 


3bb 


2m 3ab S ab> 


1 + + (m 


3cc 


+ m 


3dd 


+ 


2m„ ,S , 
3cd cd 


>( l + S ab) 


( 12 ) 
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N 2 = 4 [ m 3ac (s ac + S ab s bc + S ed s ad + S ab S cd s bd> + m 3ad< s ad + S ab s bd + S cd s ac 


+ S ab S cd s bc ) + m 3bc (s bc + S ab s ac + S cd s bd + S ab S cd s ad> 
+ m 3bd* s bd + S ab s ad + S cd s bc + S ab S cd s ap'] ^ 13) 

N 3 ’ 2 [ m 3aa( s bd + s bc + 2S cd s bc s bd) + m 3bb( s ad + s ac + 2S cd s ac s ad) + m 3cc( s bd * s ad 


+ 2E ab s ad s bd) + m 3dd(®bc + 3 ac + 2S ab s ac s bc)] * 4 { m 3ab[ s ad s bd 
+ s ac s bc + S cd< s ac s bd + s ad s bc>] + m 3cd[ s bc s bd + s ac s ad 

- S ab< s ac s bd + s bc s ad>]} < 14 ) 

E e 4 ( S ?d + ‘X 1 + *4) - 4 [ S ab (s ad s bd + s ac s bc> + S cd< s bc s bd + s ac s ad>] ' 2 ( s bd + s ac 

+ s bc + s ad) - 4S ab S cd (s bc s ad + s ac s bd> < 15 > 
3 4 

In equations (14) and (15), terms of order s and s were neglected because they were 

found to have negligible effect on u or its derivative with respect to internuclear 

ax 3 

distance for any configuration or intermolecular distance in this report. 

Following references 5 and 6, Rosen-like (ref. 13) orbitals were used, of the form 


1/2 


a(l) 


5 ab' 


-^h r 


e >ab al (l * « ab r al cos d abl + X ac r al cos 0 acl + X a(J r al cos 9^) 


(16) 


where 0 ab j is the angle at proton a between proton b and electron 1, etc. , is the 

distance from proton a to electron 1, « ab is an attractive distortion parameter (always 
positive) for an isolated Hg molecule a-b, and X &c is a repulsive distortion parameter 
(always negative) for the two repelling atoms a and c. The terms in equation (16) con- 
taining /c ab , X &c , and X ad each have absolute values much less than 1, so a(l) is a 
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slightly distorted Is atomic orbital. If a (1) terms containing k , /cX, and X are 
neglected, equation (10) is satisfied. The orbital exponent £ a ^ was assumed to have the 
same value as for an isolated Hg molecule a-b (ref. 5 or 13) with internuclear distance 

r ab- 

The distortion parameters X ac , and X a£ j were assumed to depend only on the 

distances r &b , r &c , and r a( j, respectively. The attractive distortion parameter for an 
isolated Hg molecule can be found from reference 5 or 13. The repulsive distortion 
parameter X„_, was calculated by a variational method for two repelling H atoms a and c 
by assuming orbitals of the form 



—■A 1/2 _7 


?r 


al 


& + X ac r al cos 0 acl> 


(17) 


The required variational method and integrals were given by Rosen (ref. 13) for molecu- 
lar Hg, so for two repelling H atoms it was necessary to change the signs of all Rosen's 
exchange terms and correct typographical errors. In Rosen's method both £ and X 


ac 


were optimized. In our version this would cause ^ * ? ab' To avoid this in calculating 

X„„ for use with H 0 -H 0 , £ was not optimized but was fixed at the value for the isolated 
cLC z z 

Hg molecule. 

Integrals of the type S . and s are known as overlap integrals. Because even- 

ao ac 

tually equations (11) to (15) had to be differentiated with respect to r ^ while holding 

r , constant, the case C , * if , had to be considered for s„ but did not affect S , . 
cd ^ab ^cd ac ab 

Hence, types S.- and s„„ are treated separately below. However, in both types, 
ao ac 

terms containing kX, k%, x 2, and X were neglected compared to terms containing k. 
Using spheroidal coordinates (ref. 14), integrals of type S ab were evaluated. 


S ab = ( 1 + K 


ab r ab 




(18) 


where 


= if r 
ab s ab ab 


(19) 


Prolate spheroidal coordinates were also used for integrals of type 
law of cosines from spherical trigonometry was also required. Let 


s ac > but the 
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°ac g r ac^ab + ^cd^ 


( 20 ) 



(21) 


t = — r (if -if ) 
ac g ac'^ab ^cd' 


Then 


ac 


- \ «ab?c/ /2 4 e ’ < ’ aC ( 2M ( c 'ac) + r ac{"ab cos ( >abc[ M < c ’ac> ' ^ac- T ac>] 

+ <t cd cos e cda [M(cr ac ) + L(o ac , T J 


where 


( 22 ) 


m(ct) s J_ + JL + JL 

3 2 3 a 


(23) 


L(c r,r ) = (— + — + + -i-|r 


5ct 


15a/ 


(24) 


and where, after integration, exponentials with arguments ±r were expanded in power 

o clC 

series, and then terms of order r* and higher were neglected. Equations (1) to (24) 

dv 

are equivalent to those in reference 6. 

For the integral m 3aa (eq. (7)), polar coordinates with origin at proton a and polar 
axis parallel to the x^ axis were used. Terms containing k?, kX, and X^ were ne- 
glected. It was necessary to use the law of cosines from spherical trigonometry to get 


m 


3aa 


2 

^ab 


Kb cos 0 abx 3 + A ac cos e acx 3 + X ad cos 0 adx 


(25) 


where is the angle at proton a between proton b and a line through proton a par- 

allel to the x g axis. Equations for m 3bb , m 3cc , and m 3dd were deduced from equa- 
tion (25). 

For integrals of the type m^^ (eq. (8)), prolate spheroidal coordinates and the law 

of cosines from spherical trigonometry were necessary. Terms involving k?. kX, and 
2 

X were neglected. The relation 


_ x 3al + x 3bl 
x 3abl ~ 


(26) 
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was needed. The result was 


1 „3 5 

m 3ab 4 ^ab r ab e 


^facf** 5 6 abc cos 6 ^ U Kb1 

+ sin 0 abc sin « abX3 cos <P cabx ,T(,7 ab )] + * ad [cos 0 abd cos e^Ufo^) 

+ sln "abd sin «abx 3 cos fdabx 3 T ( a ab>] + x bc[- cos Vc' coe ' “abx^'W 
+ sin e bac sine abX3 cos «> cabX; T (" ab >] * x bd [‘ cos "had cos ‘W^ab* 


+ sln Vd sin "abXj cos « , dabX3 T<< ’ab>]} 


(27) 


where (P ca ^ x is the dihedral angle between plane cab and a plane containing the line 


ab and the axis and where 


T(a) = — + — + —5- + 4 


5 4 ,3 1( -2 

a a 5a 1 5a 


(28) 


U(a) =i + i + — + + — 


a 5 <r 4 5a 3 15 ct 2 15gr 


(29) 


Integrals of the type rrig^ (eq. (9)) were treated like type except for the 
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orbital exponents, which may be unequal. The result was 


m 


cos 6 L (o , t ) + X 

acx„ v ac’ ac'l _ ae| 

6 vr ac 


= 1, (r p \ 3/ 2 5 acp-2 
Jac 4 paired/ r ac e l 

+ K ab{ cos S abc 005 "ac^Kac* ' L(<T aC T ac>] + Bin6 abc Bin e acx 3 cos 'Vcx., T <‘ T ac>j 
+ K od{ cos 9 cad COB e acx 3 [- U < a ac> ’ ^ac" T ac>] + Bin<, cad slne acx 3 COB “’daox^ac} 
+ \,d{ cos e acd COB 9 acx 3 [ u ("ac> ‘ L < 0 ac’ T ac>] + stne acd sl,1<, acx 3 c ° s ^daox^ac*} 
+ A bc{ cos «cab COB e acx 3 [- u <' 7 ac> " L(a ac T ac>] 


+ sin 0 cab sin 0. lcX3 cos <f> bacX; T (<7 ac ) 


I 


(30) 


where, after integration, exponentials with arguments ±t were expanded in power 

O CL 

series and then terms of order t^ c and higher were neglected. The angle <p^ acx is 

the dihedral angle between plane bac and a plane containing line ac and parallel to the 
Xg axis. 

Overlap contribution to deri vative of dipole moment . - In calculating the pressure- 
induced vibrational absorption, the partial derivative of the Xg component of the dipole 
moment with respect to the internuclear distance r . is the most important variable. 

If we let a prime indicate 3/dr from equation (11) 


N’i + Ng + Ng (Nj + N 2 + Ng)E* 

= " 


E 


(31) 


E 


The terms Nj, Ng, Ng, and E' were found by differentiating equations (12) to (15), but 
the results are too lengthy to reproduce here. To find Nj, Ng, Ng, and E r , the deriv- 
atives of the molecular integrals were required. These were found analytically from 
equations (18) to (25) and (27) to (30) and are also too lengthy to reproduce here. In 
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general, in taking the derivative with respect to r &b it should be noted that £ ab k^. 


r ac’ r ad’ 


be’ 


r bd' X ac’ X ad’ X bc’ X bd’ 


and most of the 0's and cp's with letter sub- 
scripts are functions of r^. 

A comparison of the assumptions made in references 5 and 6 and in this report in 
calculating n' is presented in table I. 


TABLE I. - COMPARISON OF ASSUMPTIONS MADE AND CONFIGURATIONS USED BY 


VARIOUS INVESTIGATORS IN CALCULATING OVERLAP CONTRIBUTION TO 
DERIVATIVE OF DIPOLE MOMENT IN H 2 - Hg COLLISIONS 



Van Kranendonk 
and Bird 
(ref. 5) 

Britton and Crawford 
(ref. 6) 

This report 

Antisymmetric system wave function, ^ 

No 

Yes 

Yes 

Only Xg component of p^ 

Yes 

Yes 

Yes 

Type mg ac integrals, as well as type 

No 

No 

Yes 

m 3ac’ retained in N 2 




NJj included in p* 
i ax 3 

No 

No 

Yes 

Terms of order s 2 retained in final 

No 

No 

Yes 

expression for E 




(IN^ + N 2 + N 3 )E t /E 2 retained in p^ 

— , 3 

No 

No 

Yes 

Value of orbital exponent £ in calcula- 

1.000 

1.000 

1. 174 

lation of X, bohr" 1 (m“*) 

(1. 890xl0 10 ) 

(1.890X10 10 ) 

(2.219X10 10 ) 

Factor of 1 + xr in integrals of type S ab 

No 

Yes 

Yes 

Complete expression for integrals of type 

No 

No 

Yes 

m 3ac and m 3ac’ exce P* for neglect 
of terms of order t 2 c and higher 




r ac r ad’ r bc and r bd included in 

No 

No 

Yes 

derivatives of integrals 




X T included in derivatives of integrals 

No 

No 

Yes 

9 T and cp T (with letter subscripts) in- 
cluded in derivatives of integrals 

No 

No 

Yes 

Number of planar configurations included 

4 

4 

4 

Number of nonplanar configurations included 

1 

0 

1 

Equilibrium internuclear distance of H 0 

a l. 416 

a l. 416 

b l. 401446 

z 

molecules, bohr(m) 

(0.7493X10" 10 ) 

(0. 7493x10 10 ) 

(0. 741-599'xl0‘ 10 ) 

°2-222 and D 222-2 included in ^ 

0 

expansion 

No 

No 

Yes 


Equilibrium internuclear distance for Rosen model (ref. 13). 
Equilibrium internuclear distance from spectroscopic data (ref. 16). 
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Configurations and expansion coefficients . - In calculating the pressure- induced 
vibrational absorption, n' is always expanded in spherical harmonics Y, (refs. 5 

Xn tU 

and 6). 6 

^rZ lit (32> 

^l - 0 l 2~° l J ‘2~~ l 2 

where 


Y lfl (6, <p) = (2ir) 


_1/2 0 7 Me 1 * 4 * 




(33) 


and 0^(0) are normalized associated Legendre functions tabulated by Pauling and 
Wilson (ref. 12) for ^2 0. For // < 0, we take 0^ = (-1)^ 0^ | |. The fact that 

molecules a-b and c-d were homonuclear required that many of the D's be 0 because of 


symmetry. In addition, we restricted our calculations to the five configurations shown 
in figure 2. Thus, the only nonzero D's that could be determined were D 0 000’ D 2000’ 


Side view 


End view 



X 1 

x 2 * O 


(a) Linear configuration 1. 


T 

> O 

(> J 

) 6 


O O 


(b) Quadrilateral configuration 2. 


O 

A 


(c) Planar perpendicular configuration 3. 
O 

O O 

O 


o 

() 

o 


9 


6 


(d) Planar perpendicular configuration 4. 


O 



(e) Nonplanar perpendicular configuration 5. 

Figure 2. -Configurations used in calculations for H 2 -H 0 collisions. 
Intermodular distance R had values from 4 to 8 bohr (2. 117xl0“ 10 
to 4. 233x10" 10 ml. The molecules all had internuclear distances of 
L 401446 bohr (OL741599xlO' 10 m), which is the equilibrium value (ref. 161. 
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TABLE H. - POLAR AND AZIMUTHAL ANGLES FOR FIVE 


CONFIGURATIONS SHOWN IN FIGURE 2 


Configuration 

Polar angle of 

Polar angle of 

Azimuthal angle 

Azimuthal angle 


molecule 

molecule 

of molecule 

of molecule 


a-b. 

c-d, 

a- 

-b. 

c- 

d, 


e r 

e 2 9 

<p v 

4*2’ 


deg 

deg 

deg 

deg 

1 

0 

0 

0 

0 

2 

90 

90 





3 

0 

90 





4 

90 

0 





5 

90 

90 



90 


D 0020 ’ D 2020 ’ D 2-222 
the configuration. 


and £> 222-2 two are ec l ua ^)- Let J be 

Then from equations (32) and (33) and table II, 



for 



(v-'x 




" D 0000 + ^ D 2000 ^ 


5 

+ '» / “ JJ 0020 + g D 2020 


\ D 0000 D 2000 ~Y^ D 0020 + | D 2020 


2 D 0000 + y~ D 2000 " fi, D 0020 " \ D 2020 

1 IT /s 5 

2 D 0000 " V— D 2000 + y 4 °0020 " ^ D 2020 

1-0000 -^000 -^ 


5 n 

0020 + g 2020 


+ — D 


2-222 


> 


15 

8 


D 2-222 


where the equality of D 2-222 anc * ^222-2 was invoIse d. Equations (34) were solved 
for D 


(34) 
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In calculating pressure- induced translational or rotational absorption with any 

model, or pressure- induced vibrational absorption with anharmonicity and H 2 

vibration- rotation interaction, » is also expanded in spherical harmonics 

x 0 


M, 


= Z L 2 Z 


) (36) 


1 1=0 


1 Z 2 '° ^ 2"~ l 2 


Due to the symmetry required by homonuclear diatomic molecules and due to the limita- 
tion to the five configurations in figure 2, the only C's that might be nonzero that can 
be determined are Cqqqq, C 20 q 0 , Cq 020 , C 202 q, C 2 _ 222 , and C 222 2 (the last two 

are equal). By symmetry, n is zero for configurations 1, 2, and 5 and 

x 3 



(37) 


Combining equations (33), (36), and (37), using table II, and invoking the equality of 
C 2-222 and C 222-2 gave 
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0 - \ C 0000 + C 2000 + y^ C I 




+ 5 C 

0020 + 2 '"2020 


' + — C , 15 p 

0020 + g *"2020 + g *" 2-222 


) 3 = 2 C ° 000 + |4 C 


5 

C 0020 " “ C 2020 


■ ( m 4 


) 3 ~ 2 c ° oo ° C2oo ° + ^4 C| 


5 

'0020 ‘ ^ C 2020 


15 


0 " \ C 0000 " C 2000 " Y~ C 0020 + “ C 2020 ~ ~ C 2-222 


( 38 ) 


analogous to equation (34). The solution of equations (38) is Cqqqq = C 20 2o = ^2-222 = 0 
and 


'2000 


3V5 


li. 


x, 


= -C 


0020 


(39) 


Although equations (32) to (39) were derived for nL and n , they are linear and 

x 3 x 3 

therefore each is applicable to the overlap and quadrupole- induced contributions sep- 
arately, where 


^x 3 ^ax 3 + ^qx 3 

(40) 

^x 3 ~ ^axg + ^qx 3 

(41) 

: - q + c 

l l^l l 2^2 aZ lMi Z 2^ i 2 q ^lMi ? 2^2 

(42) 

^lM^2^2 Da ^1^1^2^2 + DqZ 1^1^2 M 2 

(43) 


with subscript a standing for overlap contribution and subscript q for quadrupole- 
induced contribution. 
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Quadrupole- induced contribution to dipole moment. - An H atom does not have a 


quadrupole moment. However, an Hg molecule has a permanent quadrupole moment 
because of the correlation of the electrons. This correlation was neglected in estimating 
the distortion of the orbitals in the overlap calculations. Therefore these calculations 
include neither the dipole moment induced in molecule c-d by the quadrupole moment of 
molecule a-b nor the dipole moment induced in molecule a-b by the quadrupole moment 
of molecule c-d. These quadrupole-induced contributions must thus be added to the 

overlap contributions to n (eq. (40)). The expansion coefficients for these 

x 3 

quadrupole-induced contributions were given by Colpa and Ketelaar (ref. 15), and for the 
Xg component are (after correcting a typographical error and allowing for the difference 
in definition of 0^) 


C q0000 " 0 


'q2-121 7 ( Q ab A cd ' Q cd A ab^ 

5R* 


'q2000 


_ 6\/5 Q ab a cd 
5 R 4 


^q212-l a ^ab A cd ” ^cd A ab^ 
5R* 


> (44) 


'q0020 


6\/5 Q cd a ab 


5 R 4 


C q2-222 “ 0 


C q2020 7 ( Q ab A cd " Q cd A ab^ C q222-2 “ 0 

5R* 


where is the scalar quadrupole moment of molecule a-b defined by 

Q ab = “ Q XX = -Q YY = g Qzz ^ 


Here X, Y, and Z are Cartesian coordinates with origin at the midpoint of the line 
connecting a and b, with Z running along line ab. The quantities Qyy, and Q zz 

are elements of the quadrupole moment tensor of molecule a-b. The quantities a and 
A ab are avera S e polarizability and anisotropy of the polarizability, respectively, of 
molecule a-b. Equations (44) may all be written 
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c _ C * 1 ^ 1*2 ^2 

qZlM^ 2 M2 r 4 


( 46 ) 


where each c, , is a constant. 

Quadrupole- induce d contribution to derivative of dipole moment. - The expansion 
coefficients for the quadrupole- induced contribution to the Xg component of the dipole 
moment were found by differ entiating equations (44) with respect to r afe . For r ^ = r fid 
they are 


D q0000 " 0 


_ 6\^5 Q-a 


q2000 5 4 

R 


D 


q0020 


_ 6\/5 Qa?* 


5 R 4 


D 


_ 4 Q’A - QA’ 
q2020 5 4 

R 


r. _ 2 Q’A - QA’ 

12-121 -j r4 


t\ _ 2 Q’A — QA’ 

q212-l _ 7 4 

D R 


D q2-222 ~ 0 


D q222-2 “ 0 


(47) 


The nonzero expressions for DqQOOO * n reference 6 are incorrect. Equations (47) may 
all be written 


D = *1 M 1* 2M 2 

q*lM 1*2^2 r 4 


(48) 


where each d, , is a constant. 
1 l^l l 2^2 


Results and Discussion 

Repulsive distortion parameter . - This parameter was calculated with a digital 
computer as discussed in the section Analysis. The orbital exponent was fixed at 
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1. 174 bohr"* (2. 219x10*® m“*), which is the value for an H 0 molecule with internuclear 

- 1 0 " 

distance of 1. 401446 bohr (0. 741599x10 m) according to Rosen’s model (ref. 5). This 
internuclear distance is the Hg equilibrium internuclear distance obtained from spectro- 
scopic data (ref. 16) and was used for consistency with the integrated absorption coeffi- 
cient calculations (see comparison with experiment). The results for the repulsive dis- 
tortion parameter are given in table in. 


TABLE in. - REPULSIVE DISTORTION PARAMETER 
FOR TWO REPELLING H ATOMS 


Internuclear 

Orbital exponent of both atoms, £ 

distance, 

r 

1.087 bohr -1 

2. 054x10*° m' 1 

1. 174 bohr' 1 

2. 219X10 10 m" 1 



Repulsive distortion parameter, 

X 



Used in H^-H calculations 

Used in Hg-H 

2 calculations 

bohr 

m 

bohr” 1 

-1 

m 

bohr' 1 

m' 1 

2. 5 

0. 1323X10 y 

, 

-0. 8179X10 1 

-0. 1546x10*° 

-0. 9688x10' 1 

-0. 1831xl0 10 

3.0 

. 1588 

-.4760 

-. 8995X10 9 

-. 5639 

1066 

3. 5 

. 1852 

-. 2716 

5133 

-. 3170 

-. 5991x10° 

4.0 

. 2117 

1509 

-.2852 

-.1710 

-. 3231 

4. 5 

. 2381 

8114x10” 2 

1533 

-. 8819X10' 2 

-. 1667 

5. 0 

. 2646 

-.4217 

-. 7969X10 8 

-. 4352 

-.8224X10 8 

5. 5 

. 2910 

-.2120 

4006 

-.2063 

-. 3899 

6.0 

.3175 

1034 

1954 

-. 9440x10° 

-. 1.784 

6. 5 

. 3440 

-.4905x10'° 

-.9269X10 7 

-.4190 

-.7918X10 7 

7.0 

. 3704 

-.2273 

-. 4295 

-. 18 12 

-. 3424 

7. 5 

. 3969 

-.1032 

-.1950 

-.7668X10" 4 

-. 1449 

8. 0 

.4233 

-. 46 02x1 O' 4 

-. 8697x10° 

-. 3181 

-.6011X10 6 

8. 5 

. 4498 

-.2020 

-. 3817 

-.1298 

-. 2453 

9.0 

. 4763 

-.8747X10' 5 

-.1653 

-. 5206x10" 5 

-.9838X10 5 

9. 5 

. 5027 

-.3744 

-. 7075x10° 

-,2074 

-. 3919 


Overlap contributions to n and n’ . - These contributions were calculated 

x 3 21 

from equations (11) to (15), (18) to (25), and (27) to (31) by using a digital computer. 

The Hg internuclear distance had the value of 1. 401446 bohr (0. 741599x10"*® m) for the 
reasons previously given. The values of k and £ were calculated from reference 5 
for this internuclear distance and were 0. 1205 bohr"* (0. 2277x10*® m~*) and 1. 174 
bohr"* (2. 219x10*® m~*), respectively. The values of k ' and were assumed to be 
the same as for the internuclear distance of 1. 416 bohr (0. 7493x10"*® m) and were 
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TABLE IV. - EXPANSION COEFFICIENTS FOR x 3 COMPONENT OF OVERLAP DIPOLE MOMENT AND ITS DERIVATIVE IN H 2 -H 2 COLLISIONS 


Intermolecular 

distance, 

R 

Expansion coefficient for 
overlap dipole moment, 

°a2000 

Expansion coefficients for derivative of overlap dipole moment 

D a0000 ] 

D a2000 

1 

D a0020 

D a202G 

D a2-222 and D a222-2 

bohr 

m 

au 

C m 

au 

C 

au 

C 

au 

C 

au 

c 

au 

c 

4.0 

0 . 2117X10' 9 

- 0 . 5942x10" 3 

- 0 . 5038xl0" 32 

0.1450 

0 . 2323X10" 19 

0.0120 

0 . 0192X10" 19 

0.0146 

0 . 0234X10" 19 

0.0015 

0 . 0024x10" 19 

0.0002 

0 . 0003X10" 1 

4.2 

.2223 

-.4791 

-.4062 

.1139 

.1825 

.0096 

.0154 

.0119 

.0191 

.0012 

.0019 

. 0001 

.0002 

4.4 

. 2328 

-. 3768 

-.3194 

. 8895X10" 1 

.1425 

I . 0767x10" 1 

.0123 

. 0970x10 * 

.0155 

. 0099x10" 1 

.0016 

. 0009x10" 1 

.0001 

4.6 

.2434 

-.2928 

-.2482 

,6896 

.1105 

.0611 

.0098 

.0783 

.0125 

.0080 

.0013 

.0006 

.0001 

4.8 

.2540 

-.2255 

-.1912 

.5310 

. 8507X10" 20 

.0484 

. 0775xl0" 20 

.0628 

. 1006X10" 20 

.0066 

. 0106X10" 20 

.0004 

. 0007x1 0 _; 

5.0 

.2646 

-.1729 

-.1466 

.4060 

.6505 

.0380 

.0609 

.0499 

.0799 

.0054 

.0087 

.0003 

.0005 

5.2 

.2752 

-.1323 

-.1122 

.3083 

.4939 

.0296 

.0474 

.0394 

.0631 

.0044 

.0070 

.0002 

.0004 

5.4 

.2858 

-.1013 

-. 8588X10" 33 

.2325 

.3725 

.0229 

.0367 

.0308 

.0493 

.0036 

.0058 

.0002 

.0002 

5.6 

. 2963 

-. 7738X10" 4 

-.6560 

. 1742 

.2791 

.0175 

.0280 

1 .0239 

.0383 

.0029 

.0046 

.0001 

.0002 

5. 8 

.3069 

-. 5900 

-. 5002 

.1298 

.2080 

.0133 

.0213 

.0184 

.0295 

.0023 

.0037 

.0001 

.0001 

6.0 

.3175 

-.4505 

-.3819 

. 96 08x10" 2 

.1539 

. 0997X10" 2 

.0160 

. 1404x10“ 2 

.0225 

. 01 84x1 O' 2 

.0029 

. 0005X1 0" 2 

.0001 

6.2 

.3281 

-.3421 

-.2900 

.7075 

. 1133 

.0743 

.0119 

.1064 

.0170 

.0144 

• 0023 

.0003 

.0001 

6.4 

.3387 

-.2589 

-.2195 

. 5182 

. 8302X10" 21 

.0550 

. 0881X10" 21 

.0800 

. 1282X10" 21 

.0112 

. 0179X10" 21 

. 0002 

. 0004x10 " 

6.6 

. 3493 

-.1956 

-. 1658 

.3776 

.6050 

.0403 

.0646 

.0598 

.0958 

.0087 

.0139 

.0001 

.0002 

6.8 

.3598 

-.1467 

-.1244 

.2739 

.4388 

.0293 

.0469 

.0444 

.0711 

.0066 

.0106 

.0001 

.0002 

7.0 

. 3704 

-. 10 97 

-. 9300X10" 34 .1978 

.3169 

.0212 

.0340 

.0327 

.0524 

.0050 

.0080 

.0001 

.0001 

7.2 

,3810 

-. 8143X10" 5 

-.6903 

. 1422 

.2278 

.0152 

.0244 

.0240 

.0385 

.0037 

.0059 

.0000 

.0001 

7.4 

.3916 

-.6013 

-. 5098 

. 1019 

.1633 

.0109 

.0175 

.0175 

.0280 

.0028 

.0045 

.0000 

.0000 

7.6 

.4022 

-.4426 

-.3752 

. 7273X10" 3 

. 1165 

. 0774X10" 3 

.0124 

. 1271X10" 3 

.0204 

. 0207X10“ 3 

.0033 

. 0002X10" 3 

. 0000 

7.8 

.4128 

-.3226 

-.2735 

.5172 

. 82 86x10" 2 2 

. 0547 

. 0876X10" 22 

.0918 

. 147 1X1 0' 22 

.0152 

. 0244x1 O' 22 

.0001 

. 0002 X 10" 

8.0 

.4233 

- . 2347 

-.1990 

, . 3666 

.5873 

.0384 

.0615 

, .0661 

.1059 

.0111 

.0178 

.0001 

.0001 


a ^a0020 = ” C a2000 



-0. 037 bohr -2 (-0. 132xl0 20 m" 2 ) and -0. 247 bohr” 2 (-0. 882xl0 20 m 2 ), respectively 
(ref. 5). The calculations were carried out for the five configurations in figure 2 and 


for intermolecular distances from 4 to 8 bohr (2. 117x10 
pans ion coefficients are tabulated in table IV. 


-10 


to 4. 233x10 m). The ex 


Quadrupole- induced contributions to n and p' . - These contributions required 


values for Q, a, and A and their derivatives Q', a', and A' for Hg. The most re- 
liable values for Q are the theoretical values of Kolos and Wolniewicz (ref. 9), who 
tabulated Q for a number of internuclear distances. From reference 9, Q was found 
to be 0. 45822 atomic units (au) (2. 0556xl0” 4 ^ C m 2 ), and Q' was found to be 0. 5314 au 
(4.505x10 Cm). The most reliable values for a and A are the theoretical values 
of Kolos and Wolniewicz (ref. 10), who again tabulated values for a number of inter- 
nuclear distances. From reference 10, a was found to be 5. 1849 au (8. 5482X10" 44 
C 2 m 2 J”*), a ' was found to be 4. 351 au (1. 356xl0” 3 ^ C 2 m J“ 4 ), A was 1. 8076 au 
(2. 9801X10” 41 C 2 m^J" 1 ), and A’ was 3. 382 au (1. 054x10” 30 C 2 m J" 1 ). All these 
values are for an internuclear distance of 1. 401446 bohr (0.741599x10”^ m). 

Using these values, the coefficients for the quadrupole- induced contributions to 

x 3 


and u . ' were calculated from equations (44) and (46) to (48) and are tabulated in 
x 3 

table V. The quadrupole- induced contributions to n and p' for the five configura- 


TABLE V. - EXPANSION COEFFICIENTS FOR QUADRUPOLE-INDUCED 
CONTRIBUTION TO COMPONENTS OF DIPOLE 
MOMENT AND THEIR DERIVATIVES 


Collision 

Coefficient for quadrupole- induced 

Coefficient for quadrupole- induced 


contribution to component 

contribution to derivative of com- 


of dipole moment 

ponent of dipole moment 


Coefficient 

Value of^co efficient 

Coefficient 

Value of coefficient 



au 

C m 5 


au 

C m 4 

V H 2 

c 2000 

6. 3749 

0 . 42377X10' 69 

d 0000 

0 

0 


c 0020 

- 6.3749 

-.42377 

d 2000 

7.393 

. 9287X10* 59 





d 0020 

- 5 . 350 

-.6721 

5920X10' 60 





d 2020 

-.4713 





d 2-222 

0 

0 





d 222-2 

0 

0 

h 2 -h 

c 00 

c 20 

0 

3.9124 

0 

. 26007X10* 69 

d 00 

d 20 

0 

4 . 537 

0 

. 5699X10" 59 


c 40 

g 21 

0 

- 2.2588 

0 

15015x10 *" 

d 40 

h 21 

0 

- 2.620 

0 

3291x10 ”" 


u 
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TABLE VI. - COEFFICIENTS FOR QUADRUPOLE-INDUCED CONTRIBUTION 
TO COMPONENTS OF DIPOLE MOMENT AND THEIR 


DERIVATIVES FOR VARIOUS CONFIGURATIONS 


Collision 

Configuration 

Component 

Coefficient for quadrupole- 

Coefficient for quadrupole- 





induced contribution to 

induced contribution to 





component of dipole 

derivative of compo- 





moment 3 * 

I nent of dipole 







moment 3 





au 

C m 5 

au 

Cm 4 

H 2- H 2 

1 

X 

3 

0 

0 

1. 108 

0 . 0139 X 10" 58 


2 



0 

0 

- 1.438 

-.0181 


3 



10.691 

. 71067 X 10" 69 

11. 846 

. 1488 


4 



- 10.691 

-.71067 

- 9 . 524 

1196 


5 



0 

0 

- 1.438 

-.0181 

h 2 -h 

1 

X 

1 

0 

0 

0 

0 


1 

X 

3 

6.1860 

. 41121 X 10" 69 

7 . 174 

. 9012 X 10' 59 


2 

X 

1 

- 3. 0930 

-.20560 

- 3 . 587 

-.4506 


2 

X 

3 

1 . 5465 

. 10280 

1.793 

.2252 


3 

X 

1 

0 

0 

0 

0 


3 

X 

3 

- 3. 0930 

-. 20560 X 10" 69 

- 3 . 587 

4506 x 10' 59 


El 4 

Coefficients must be divided by R to obtain the component of the dipole moment 
or its derivative. 


tions were then calculated from equations (38) and (34), respectively. These are tabu- 

4 

lated in table VI (the values given must be divided by R ). 

Total a and j±' . - These totals were calculated from equations (40) and (41) 
x 3 x 3 

for the five configurations and are plotted in figures 3 and 4. In general, they are 

strongly decreasing functions of intermolecular distance. 

Comparison with other calculations. - Present results for (eq. (40)) and fil 
- x 3 x 3 

(eq. (41)) are compared in figures 3 and 4 with other calculations for various configura- 
tions. In each figure, two other calculations are given: (1) an ab initio configuration- 
interaction calculation by Patch (ref. 11), and (2) a semiempirical overlap calculation 
(ref. 6 or 17) plus the same quadrupole-induced contributions as in this report. Use 
of the quadrupole-induced contributions of this report, which are based on references 9 
and 10, and on equations (34), (38), (44), and (47), prevents older, obviously inaccurate 
calculations of the quadrupole-induced contributions from confusing the comparison. In 
figures 3, 4(a), and 4(c), all three calculations are in reasonable agreement between 
4 and 5 bohr (2. 117x10”^ and 2. 646x10” m). The disagreement between the present 
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X 3 component of dipole moment, p x , C m 


10-29 , 


10-30 


C- =rio' 


This report (eq. (40)) 

Patch (ref. 11) 

Overlap contribution from 

ref. 17, quadrupole-induced 
contribution calculated from 
refs. 9 and 10 and from 
eqs. (38) and (44) 


10 


-31 


E 

CD 

a 

TZI 

o 

Ti 

CD 

c 

a 


£ 10 


-2 





10-32 


10 


-3 1 


4 5 6 

Intermolecular distance, R, bohr 


1.5 


I 

2.5 


1 

3.0 


2.0 2.5 3.0 3.5 

Intermolecular distance, R, m 


4. 0x10 


,-10 


Figure 3. - X 3 component of dipole moment for H 2 -H 2 configuration 3. Dipole 
moment of configuration 4 is the same except for sign. Configurations 1, 2, 
and 5 have no x 3 component of dipole moment. 
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(a) Linear configuration 1. 



Intermodular distance, R, bohr 




1.5 2.0 2.5 3.0 3.5 4.0xl0 ' 10 

Intermodular distance, R, m 

(b) Quadrilateral configuration 2 and nonplanar perpendicular con- 
figuration 5. The difference in the derivatives for configuration 2 
and 5 according to equation (41) is too small to show. 

Figure 4. - Derivative of X 3 component of dipole moment in H 2 "H 2 collisions. 


24 



Derivative of X 3 component of dipole moment, dp x Idr^, au 




1 \ 

\ 

\ 

\ 

\ 

\ 

\ 



lc) Planar perpendicular configuration 3. 



I I I J I I in 

1.5 2.0 2.5 3.0 3.5 4.0x10 iU 

Intermolecuiar distance, R, m 

Id) Planar perpendicular configuration 4. 

Figure 4. - Concluded. 
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results and Patch's ab initio calculation at other values of R has two causes: (1) the 
ab initio calculation has insufficient allowance for electron correlation in the wave 
function and no distortion of orbitals and, hence, is not valid for large R; and (2) the 
present results do not include configuration interaction and hence are not valid for 
small R. In figures 3, 4(a), and 4(c) the present results agree substantially with over- 
lap calculations from reference 6 or 17 at all values of R. 

In figures 4(b) and (d) the agreement of the three calculations is poor for all R 

values. The reason is unknown. Fortunately, the magnitudes of »' in figures 4(b) 

x 3 

and (d) are only about one-half those in figures 4(a) and (c), so the lack of agreement be- 
tween the present work and reference 11 is not too serious. The probable error should 
be reduced by fairing the results together as discussed in the next section. In addition, 
figures 4(b) and (d) show a disagreement between the present overlap contribution and 
the overlap contribution from reference 6. This disagreement is principally due to the 
different values of £ used in calculating X: 1. 000 bohr - * (1. 890x10*** m - *) in refer- 
ence 6, and 1. 174 bohr - * (2. 219x10*** m - *) in the present work. The latter value is the 
same as used for the orbitals in the Hg-Hg calculation and is thus consistent. Also, it 
gives better agreement with the measured integrated absorption coefficient of the funda- 
mental as outlined in the following section. 

Comparison with experiment. - The integrated absorption coefficient of a pressure- 
induced vibrational transition may be said to vary roughly as the square of the derivative 
of the dipole moment averaged over intermolecular distance. To provide a comparison 
of the results in tables IV and V with experiment (ref. 18), the integrated absorption 
coefficient of the Hg-Hg pressure- induced fundamental vibrational transition was cal- 
culated for a temperature of 298 K. The method of calculation was based on reference 6 
with the following major changes: 

(1) The Ho-Hr, average interaction energy was assumed to be a Morse potential 

^ ^ - 10 

(ref. 19) for R greater than 4 bohr (2. 117x10 m). For R less than 4 bohr, an ex- 

ponential was fitted to the Morse potential at 4 bohr and an ab initio energy (ref. 11) at 
2. 5 bohr (1.323x10'*° m). 

(2) ©2-222 an< * ^222-2 terms were included. 

(3) Anharmonicity and vibration-rotation interaction were included, thereby re- 


quiring c 2000 . 

(4) The quantities 


from reference 11 were used for R between 


2. 5 and 4 bohr (1. 323x10" 10 and 2. 117x10" 10 m). The quantities n and u' from 

x 3 x 3 

tables IV and V were used for R between 6 and 8 bohr (3. 175x10" 10 and 4. 233x10" 10m). 


Between 4 and 6 bohr the expansion coefficients for n x and /j' x were faired into re- 

^ 3 3 

suits from reference 11 at R = 4 bohr (2. 117x10" 10 m) and tables IV and V at R = 6 bohr 
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jf 


(3. 175x10" ^ m) by setting each expansion coefficient equal to f^R - ^ + f 2 + fgR + f^R^, 

where f 1 to f . are constants. This procedure was necessary for the reasons dis- 
’ _i o 

cussed in the preceding section. For R greater than 8 bohr (4.233x10 m), there is 

little contribution to the absorption coefficient at 298 K, so the results of table IV were 

extrapolated and added to results from table V. 

(5) Typographical errors were corrected. 

(6) The Xj and Xg components of p and m’ were neglected. 

Using these assumptions, the integrated absorption coefficient of the fundamental 
was calculated to be 87 percent of the experimental value for the binary integrated ab- 
sorption coefficient (ref. 18). The calculated value was low principally because the x^ 
component of p' was neglected (it was zero for the five configurations in fig. 2 but 
not for many others). 

Potential applications . - The H 2 ~H 2 results in tables IV and V, together with re- 
sults in references 11, 16, and 19 and available line shapes, are sufficient for calcula- 
tions of Hg-Hg pressure-induced vibrational, rotational, and translational absorption 
coefficients at temperatures up to 7000 K without recourse to the very dubious ex- 
trapolations and assumptions appearing in earlier work. 


DIPOLE MOMENT AND ITS DERIVATIVE IN H 2 "H COLLISIONS 

Analysis 

Just as for Hg-Hg, the dipole moment of Hg-H and its derivative consist of two 
additive contributions: an overlap contribution and a quadrupole- induced contribution. 
The analysis is similar to that for H 2 -H 2 except that (1) the repulsive distortion param- 
eter X is based on a repulsive H-H calculation using an orbital exponent £ which is 
the average of the orbital exponents for H and H 2 , (2) the x^ components of p and p' 
are included as well as the Xg components, (3) no molecular integrals are neglected, 
and (4) no series expansions in powers of r are used in evaluating molecular integrals 
because of the appreciable difference between H 2 and H orbital exponents. All equations 
are in atomic units (bohr, hartree, and electron charge). 

Coordinates . - The coordinate system is shown in figure 5. The origin of the 
Cartesian coordinates x^, x 2 , and Xg is on the intermolecular axis one-third of the 
way from molecule b-c to atom a. The relative position of the molecule is given by the 
polar angle 6 j, and the molecule is always in the x^-Xg plane. Hence, the x 2 com- 
ponents of p and p* are zero, and is zero. 
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c R 

3 


r R 1 

Figure 5 . - Coordinates in H 2 "H collisions. The molecule is b-c and the atom is 
is at a; protons are located at a, b, and c; the x 3 axis passes through the 
centers of the molecule and the atom; and <p ± are polar and azimuthal 
angles, respectively. 


Overlap contributio n to dipole moment. - For intermolecular distances of 4 bohr 
(2. 117x10" m ) or greater, it is reasonable to neglect configuration interaction, so the 
antisymmetric system wave function was assumed to be (ref. 20) 

* = V'l - ^2 (49) 


where the determinantal wave functions ^ and are given by 

V'J s (3 !)" 1/2 det[a(l)a(l),b(2)^(2),c(3)a(3)] 

► 

xp 2 = (3!)" 1/2 det[a(l)a(l),b(2)or(2),c(3)j8(3)] 

J 

The x. component of is given by 
J * 


M. 


- f dV l dV 2 dV 3 

J i=l 


ax j 


A 2 dv l 


(j = 1, 3) 


dV 2 dV 3 


(50) 


(51) 
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If equation (10) is assumed to hold, and integrals are defined as in equations (5) to (9) 
except for generalizing the subscript 3 to j, equation (51) becomes 


? 


where 

N Xj * ~ 2 ( m jaa 


Max, = 0 =1, 3) (52) 

J 

+ m )bb + m jcc + m jaa®bc * 2m jbc s bc) + 2 ( m jac s ac * m )ab s ab + m )ab s bc l, ac 

+ m Jac S bc s ab) + m jbb s ac + m )ec4 + 2m jbc s ab s ac 0 ' *< 3 > < 53 > 

E * 2 + 2 4 - 4 - 4 - 2S bc s ab s ac < S4 > 


As for Hg-Hg, Rosen-like orbitals were used: 

, 3 \l/2 

c\ . 


a(i) = — e a ai (l + X ab r ai cos 0 &bi + A ac r ai cos 0 &ci ) 


1/2 


b(i) = (^] e ?bcr>bl (l + Kr bi cos 0 bc - + A^. cos 0 feai) 


r(i= 1,2,3) ( 55 ) 


c(i) = 


- \ 1/2 

— | e ?bcrci (l + Kr ci cos 0 cbi + A ac r c . cos 9^) 

77 


where if is the orbital exponent for the atom (1 bohr" 1 or 1. 890xl0 10 m -1 ) and k is 
the attractive distortion parameter for an isolated Hg molecule with internuclear dis- 
tance r bc> 

As for Hg-Hg, there are two kinds of overlap integrals. The integral S bc was 

evaluated by equations (18) and (19). Because of the appreciable difference between 

Hg and H orbital exponents £ bc and £ , respectively, integrals of type s ac were 

derived without truncated series expansions in powers of r , where 

ac 
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( 56 ) 


T sir (C - if ) 
ac 2 ac Vb a b bc' 


The result is 


_lL K \ 3/2 r 3 e - a ac 
s ac g Y a ^ bc ) r ac 


i 2M(a . t ) + r /c cos 0 , M(cr ,r ) + L(ct , r J > 

^ v ac’ ac' ac cbaL ac’ ac' ' ac’ ac'JJ 


where 


( 57 ) 


°ac 2 r ac^a + ^bc^ 


( 58 ) 


M (a, r) h I [A 1 (r)J 3 (a) - AgfrJJ^a)] 
L(a,r) =1 [A 4 (r)J 2 (a) - A 2 (t)J 4 ((t)] 


( 59 ) 

( 60 ) 


A n ( T ) H " e " 


1 + n - 1 + (n - l)(n - 2 ) + 


, (n - 1): , (n - 1): 
n-1 n 

T T 


(-D 1 


1 ) , (-l) u ~*(n- l)(n - 2) , (n-1); (n-1)! 


x n-2 


+ • « . + 


n-1 n 

T T 


( 61 ) 




n 


a 


n- 1 


1 + a + — + 
2! 


(n - 1)! 


( 62 ) 


Integrals of type m^ aa were treated in a fashion similar to m 3aa for with 

the result 


m jaa = “ ( x ab cos 0 abx. + X ac cos 0 ac Xj ) 0 b1 - 3 > 

^a 


m jbb “ 


r 2 

^bc 


( K cos 0 bcx j - X ab cos 0 ab Xj ) 0 = X ’ 3 > 


( 63 ) 


( 64 ) 
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The integral m j cc is similar to m- bb . 

The integral m^ bc was derived just as for Hg-Hg 


m jbc 4 **bc r bc e 


I c °s e bca cos 9 V (a bc ) 

J 


derived* 

"“{‘4 

+ sin 9 bca sin 0 bcXj cos (» abcx T(o bc )l + X a< |- 

i T( ^ ) ]} 0 ' 1 ’ 


+ sin e cba sin 0 bcXj cos ? abcx 


cos « cba nos 9 bcXj U(o b( .) 
3 ) ( 65 ) 


where cr bc , T(a bc ), and U(cr bc ) were found from equations ( 19 ), ( 28 ), and ( 29 ), respec- 
tively. 

Just as for s„ , integrals of type m. were derived without series expansions in 

aC ]aC 

powers of t 


ac 


\3/2 5 -cr : 


m ’~- = | ^ a ? bc) ' r ac e a °l - 2 cos 6 »«*Ma»n> T »o)l— + x 


]ac 


acx. ' ac’ ac'i ' "ac] 
J \ r ac 


+ *{ cos e cba cos «acXj[-^ac. T ac) ' ^ac- T ac>] 

+ sin 9 cba sin 6 acXj cos ?> bacx T(o aCj T ac )J 

+ x ab{ cos 9 abc 008 ‘ , acx j [ U < a ac>’ac> ‘ T ac>] 

> T acJ 


+ sin 0 , sin 9 cos cp , T(cr 

abc acx. ^bacx. v ac 

J J 


(j = 1 , 3 ) ( 66 ) 


where 


T(( 7 ,t)h i jj^jA^T) - A 3 (t)] + J 3 (ct)[a 5 (t) - Aj(t)J + J^A^r) - A 5 (t)]} 


( 67 ) 
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U(ct,t) =± [j 5 (ct)A 3 (t) - J 3 (a)A 5 (T)] 


( 68 ) 


and L(cx , t ) is found from equation (60). 

3.C clC 

Overlap contribution to derivative of dipole moment. - In calculating the pressure- 
induced vibrational absorption, the partial derivative of ~j± with respect to the inter- 
nuclear distance r bc is necessary. If we let a prime indicate 3/ 3r bc , from equa- 
tion (52) 


N T N E' 
x. x. 

n ' = 1 * 

MaX i E e 2 


(j = 1, 3) 


(69) 


where N' and E' were found by differentiating equations (53), (54), (57), (63) to (66), 

X j 

etc. , analytically. 

Configurations and expansion coefficients . - In calculating the pressure- induced 
vibrational absorption for H 0 colliding with an atom (ref. 21), u* is expanded in nor- 

i x 3 

malized Legendre polynomials 0^(8) 


Hi, - J 

a 1=0 


(70) 


Side view 



-R— — 

(a) Linear configuration 1. 


End view 


x 2 




O 


(b) Scalene configuration 2. The molecule is at 45° to the intermolec- 
ular axis (x 3 axis). 

O 

O 

o 


(c) Isosceles configuration 3. 

Figure 6. - Configurations used in calculations for H 2 -H collisions^ 
Intermolecular distance R had values from 4 to 8 bohr (2. 117x10 1 
to 4.233xl0 -10 m). The molecule had an internuclear distance of 
1. 401446 bohr (0. 741599x10“ 10 m). 
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The fact that molecule b-c is homonuclear requires that 0 ^q be zero when l is odd. 

In addition, calculations were restricted to the three configurations shown in figure 6. 
Thus, the only nonzero D f s that could be determined were Dqq, D 2 q, and D^q. From 
equation (70) and figure 6 



Equations (71) were solved for the D's as follows: 



In calculating absorption coefficients it is convenient to take linear combinations 

of u and p . For any configuration 
X 1 x 2 

H i) H + iM x 2 )’ 

M(i) = (2) ^ (^ x i + i ^ x 2 > )’ 
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where i is the imaginary unit. For the configurations in figure 6, the (1) and (-1) 
components are equal. The quantity was expanded in spherical harmonics 




m il=-i 


i & 




( 74 ) 


Because of the symmetry and the three configurations considered, all of which have 
i± x , and n x equal to zero and two of which have zero the only nonzero 

term in equation (74) that could be determined was 


H 21 = 


VS " VS ( Mx i 


( 75 ) 


In calculating pressure- induced translational or rotational absorption with any 
model, or pressure- induced vibrational absorption with anharmonicity and vibration- 
rotation interaction, ~n is needed, as well as 7 l '. Hence, u was expanded. 

' Xn 


- L WV 


1=0 


(76) 


which gave 


C 00 “ 


C 20 “ 


C 40 ” 



\ 

+ i &( 

15 \ x 3 ^ 

15 \ X 3/ 2 

5 V 

2° L ' 

1 + 16 / > 

1 " 12 

21^' 

' 1 21^15 ( X 3* 

2 7 yf \ 0 

32 L ^ 

1 - 64 (u \ 

1 + 32 

105 V 2 ' x v 

’l 105 V 2 X3 ' 

2 105 y /2 




'N 


M-v 


> 


(77) 


Also, ijl^ was expanded 


CO l 


“(D'Z X, V^G (R)Y w,,,,,) 

1=0 v=-l 


( 78 ) 
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The only nonzero term in equation (78) that could be determined from the three configu- 
rations was 


G 


21 " 


Vl5 




(79) 


Although equations (70) to (79) were derived for q' and q, they are linear and 
therefore equally applicable to the overlap and quadrupole- induced contributions to the 
expansion coefficients and components of q' and ~jj. just as for Hg-Hg. 

Quadrupole- induced contribution to dipole moment . - The quadrupole- induced con- 
tribution to the dipole moment when a homonuclear diatomic molecule collides with an 
atom can be found from Colpa and Ketelaar's (ref. 15) equations for the electric field 
strength F due to a point quadrupole and the relation u = F a , where a is the 
polarizability of an H atom. The expansion coefficients are 



(80) 


Equations (80) may be written 




, _ 10 
'qZO " "7 

R 


J q2 1 


g 21 


R' 


(81) 


where each c^q and g 21 


are constants. 
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Quadrupole- induced co ntribution to derivative of dipole moment. - The expansion 
coefficients for these derivatives are found by differentiating equation (80) 



Equations (82) may be written 



H q21 




(82) 


(83) 


where each d^Q and h,^ 


are constants. 


Results and Discussion 


Repulsive distortion para meter. - This parameter was calculated just as for H 2 ~H 2 

except that £ was set equal to the average of the H 0 and H values, namely 1. 087 bohr~l 
10-1 ^ 

(2. 054x10 m ). The resulting values of X are given in table m. 

Overlap contributions to u and ul. • - These contributions were calculated from 


-x and . 
i x i 


equations (18) and (19), (28) and (29), (52) to (54), and (56) to (69) by using a digital computer. 
Values of the H 2 equilibrium internuclear distance and k, k', , and j^ c were the 

same as used for Hg-Hg, but, of course, was taken to be 1. 000 bohr -1 (1. 890x10-^ 
m - ■*■) . The calculations were carried out for the three configurations in figure 6 and for 
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TABLE VH. - EXPANSION COEFFICIENTS FOR COMPONENTS OF OVERLAP DIPOLE MOMENT IN H 2 -H COLLISIONS 
Intermolecular Expansion coefficients for components of overlap dipole moment 


distance, ’ ; 

R C a00 C a20 C a40 G a21 


bohr 

m au 

C m 

— 

au 

C m 

au 

C m 

au 

C m 

4.0 

0.2117 xl 0“ 9 -0. 4086 x 10" 1 

-0. 3464 X 10' 30 

-0. 03 53 x 1 0' 1 

-0. 0299 X 10" 30 

- O . OOOlxlO -1 

-0. OOOlxlO ' 30 

0. 1580 xl 0“ 2 

0. 1339 X 10” 31 

4.2 

.2223 -.3285 

-.2785 

-.0292 

-.0248 

-.0001 

-.0001 

.1186 

.1005 

4.4 

.2328 -.2623 

-.2224 

-.0239 

-.0203 

-.0000 

-.0000 

. 8906 x 10“ 3 

.7550 X 10' 32 

4.6 

.2434 -.2080 

-.1763 

-.0195 

-.0165 

-.0000 

-.0000 

.6667 

. 5652 

4.8 

.2540 -.1639 

-.1390 

-.0158 

-.0134 

-.0000 

-.0000 

.4981 

.4223 

5.0 

.2646 -.1282 

-.1087 

-.0127 

-.0180 

-.0000 

-.0000 

.3708 

.3144 

5.2 

.2752 -. 9972 xl 0’ 2 

-. 8454 x 10* 31 

1007 X 10 -2 

-.0854 X 10' 31 

-.0002 x 10“ 2 

-.0002 x 10“ 31 

.2750 

.2331 

5.4 

.2858 -.7707 

-.6534 

-.0795 

-.0674 

-.0002 

-.0001 

.2032 

.1723 

5.6 

.2963 -.5924 

-.5022 

-.0623 

-.0528 

-.0001 

-.0001 

.1495 

. 1267 

5.8 

.3069 -.4527 

-.3838 

-.0484 

-.0410 

-.0001 

-.0001 

.1095 

.9283 X 10” 33 

6.0 

.3175 -.3443 

-.2919 

-.0374 

-.0317 

-.0001 

-.0001 

.7996 X 10” 4 

.6779 

6.2 

.3281 -.2606 

-.2209 

-.0287 

-.0243 

-.0000 

-.0000 

.5811 

.4926 

6.4 

.3387 -.1963 

-.1664 

-.0218 

-.0185 

-.0000 

-.0000 

.4204 

.3564 

6.6 

.3493 -.1473 

-.1249 

-.0165 

-.0140 

-.0000 

-.0000 

.3031 

.2570 

6.8 

.3598 -.1100 

-.9326 x 10" 32 

-.0124 

-. 1051 xl 0” 32 

-.0000 

-.0002 x 10” 32 

.2176 

.1845 

7.0 

.3704 -.8191 x 10” 3 

-.6944 

-.093 ixio ” 3 

-.0789 

-. OOOlxlO 3 

-.0001 

.1557 

.1320 

7.2 

.3810 -.6077 

-.5152 

-.0693 

-.0588 

-.0001 

-.0000 

.1110 

.9410 X 10" 34 

7.4 

.3916 -.4492 

-.3808 

-.0514 

-.0436 

-.0000 

-.0000 

.7890 x 10” 5 

.6689 

7.6 

.4022 -.3311 

-.2807 

-.0379 

-.0321 

-.0000 

-.0000 

.5588 

.4737 

7.8 

.4128 -.2432 

-.2062 

-.0278 

-.0236 

-.0000 

-.0000 

.3941 

.3341 

8.0 

[ .4233 -.1782 

-.1511 

-.0204 

-.0173 

.0000 

.0000 

.2779 

.2356 
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TABLE VIE. - EXPANSION COEFFICIENTS FOR DERIVATIVES OF COMPONENTS 


OF OVERLAP DIPOLE MOMENT IN Hg-H COLLISIONS 


Intermolecular 

distance, 

R 

Expansion coefficients for derivatives of components of overlap dipole moment 

D a 00 

D a 20 

D a 40 

H a 21 

bohr 

m 

au 

C 

au 

C 

au 

c 

au 

C 

4.0 

0 . 2117 x 10 ^ 

0 . 3889 X 10" 1 

0 . 6231 xl 0~ 20 

- 0 . 0157 xl 0‘ 1 

- 0 . 0252 X 10' 20 

0 . 0005 x 1 0 -1 

0 . 0008 x 10" 20 

; 0 . 1205 x 1 0" 2 

0 . 1931 X 10" 21 

4.2 

.2223 

. 3001 

.4808 

-.0149 

-.0239 

.0003 

. 0005 

. 9060 X 10 -3 

. 1452 

4.4 

.2328 

.2299 

. 3683 

-.0136 

-.0218 

.0002 

. 0003 

.6797 

.1089 

4.6 

.2434 

.1748 

.2800 

-.0121 

-.0194 

.0001 

.0002 

. 5075 

. 8130 X 10" 22 

4.8 

.2540 

.1319 

.2113 

-.0106 

-.0170 

.0000 

.0001 

.3772 

.6043 

5.0 

.2646 

. 9886 X 10" 2 

. 1584 

-. 0908 X 10' 2 

-.0145 

, . 0002 X 10" 2 

. 0000 

.2790 

.4470 

5.2 

.2752 

.7355 

.1178 

-.0767 

-.0123 

-.0001 

-.0000 

.2052 

.3288 

5.4 

.2858 

. 5433 

. 8704 xl 0“ 21 

-.0638 

-. 1022 xl 0” 21 

-.0002 

-. 0002 X 10" 21 

.1500 

.2403 

5.6 

.2963 

.3986 

.6386 

-.0524 

-.0840 

-.0002 

-.0003 

.1090 

.1746 

5.8 

.3069 

.2904 

.4652 

- . 0426 

-.0682 

-.0002 

-.0003 

. 7868 X 10" 4 

.1261 

6.0 

.3175 

.2102 

. 3368 

-.0342 

-.0548 

-.0002 

-.0003 

: . 5642 

. 9039 X 10" 23 

6.2 

.3281 

.1511 

.2421 

-.0273 

-.0437 

-.0002 

0002 

.4020 

.6440 

6.4 

. 3387 

.1079 

. 1729 

-.0215 

-.0344 

-.0001 

-.0002 

.2837 

.4545 

6.6 

.3493 

. 7657 X 10" 3 

.1227 

-. 1687 xl 0~ 3 

-.0270 

- . OOlOxlO -3 

-.0002 

.1990 

.3188 

6.8 

.3598 

.5394 

. 8642 X 10" 22 

-.1312 

-. 2102 X 10' 22 

-.0007 

r . 0011 xl 0" 22 

.1383 

.2216 

7.0 

.3704 

.3776 

.6050 

-.1013 

-.1623 

-.0005 

-.0009 

. 9537 X 10" 5 

.1528 

7.2 

.3810 

.2625 

.4206 

-.0777 

-.1245 

-.0003 

-.0005 

.6519 

. 1044 

7.4 

. 3916 

.1811 

.2901 

-.0592 

-.0948 

-.0002 

-. 0004 

.4403 

. 7054 X 10" 24 

7.6 

.4022 

.1240 

. 1987 

-.0449 

-.0719 

-.0001 

-.0002 

.2943 

.4715 

7. 8 

.4128 

. 8415 X 10" 4 

. 1348 

-. 3386 X 10' 4 

-.0542 

-. 0005 xl 0~ 4 

-.0001 

. 1945 

.3116 

8.0 

.4233 

. 5662 

. 9071 X 10" 23 

-.2541 

-. 4071 X 10' 23 

-.0002 

-. 0002 X 10" 23 

.1267 

.2030 



Component of dipole moment, p x> , C m 






intermolecular distances from 4 to 8 bohr (2. 117x10" 10 to 4. 233x10" 10 m). The expan- 
sion coefficients are tabulated in tables VII and VTIL 

Quadrupole- induced contributions to n x an£ i M x / " The same values of Q and 

Q’ were used as for Hg-Hg. The polarizability a a of the H atom was obtained from 
Pauling and Wilson (ref. 12) and had the value 4. 5 au (7. 419x10 ^ m^ J ^). 

Using these values, the coefficients for the quadrupole-induced contributions to 
n n , n' , and U* were calculated from equations (80) to (83) and are tabulated 
X 1 x 3 X 1 x 3 

in table V. The quadrupole-induced contributions to M x » M x > M x > and M x for the 


L 1 x 3 1 3 

three configurations were then calculated from equations (71), (73), (74), (76), and (78) 

and are tabulated in table VI ( the values given must be divided by R 4 ). 


Total fx. 


x. 


M-v 




' , and /i’ . - These totals were found by summing the overlap 


x 


and quadrupole-induced contributions for the three configurations and are plotted in 
figures 7 and 8. In general, their magnitudes are strongly decreasing functions of inter- 
molecular distance. 



Intermolecular distance, R, m 

Figure 8. - Derivatives of Xj and components of dipole moment in FU-H collisions. The derivative 
of the x 1 component of dipole moment is zero for linear configuration 1 and isosceles configuration 3. 
Overlap and quadrupole-induced contributions are included. 
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Comparison with other investigations . - For the range of intermolecular distance 
used, there are no published theoretical or experimental values of "Jl or ~jl' . 

Potential applications . - No application is foreseen until ab initio calculations are 
made for R less than 4 bohr (2. 117x10”^ m). The semiempirical method in this re- 
port is not applicable for R less than 4 bohr because it neglects configuration inter- 
action. When used with ab initio calculations and a realistic average interaction energy 
(ref. 22), the Hg-H results in tables V, VH, and vm should make possible the calcula- 
tion of the Hg-H pressure-induced vibrational absorption coefficient. This coefficient 
has not previously been calculated. 


CONCLUDING REMARKS 

The electric dipole moment and its partial derivative with respect to Hg internu- 
clear distance were calculated for Hg-Hg and Hg-H collisions by using a semiempirical 
theory. Four planar configurations and one nonplanar configuration were employed for 
H 0 -H o . Three configurations were used for H 0 -H. Intermolecular distances ranged 

“ * -10 _ t n “ 

from 4 to 8 bohr (2. 117x10 to 4. 233x10 m). For intermolecular distances less 

than about 4 bohr, the semiempirical theory is not valid because of neglect of configura- 
tion interaction. To overcome this obstacle, the dipole moment and its derivative for 
Hg-Hg from this report were faired into values from previous ab initio calculations for 
small intermolecular distances. The resulting integrated absorption coefficient at 
298 K for the fundamental pressure- induced vibrational transition was 87 percent of the 
experimental value. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 19, 1970, 

122-29. 
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APPENDIX -SYMBOLS 


Two sets of units are given for many symbols in this appendix: atomic units and 
SI units. In the atomic units, "charge" means "charge of an electron. " In the text, 
tables, and figures, atomic units are frequently abbreviated au. All equations in the 
main text are in atomic units. 


A n (r)(n =1,2, . . . , 5) 
a(i)(i= 1,2, 3,4) 
b(i)(i= 1,2, 3, 4) 

C aZ0 (Z = °> 2 ’ 4) 

C aZl^l^2M2 
C lQ (l = 0, 2, 4) 

C ql0 (l = 0,2,4) 

C q^lM^2^2 

c(i) (i - 1,2, 3, 4) 
c l0 (l = 0, 2, 4) 

D aZO a = 0, 2, 4) 

D aZ 1 ^ 1 Z 2 M2 


functions of r 

orbital of electron i about proton a, bohr - *^; 
orbital of electron i about proton b, bohr - '^; 
expansion coefficient for a in H 0 -H collision, charge 

oXg Z 

bohr; C m 

expansion coefficient for in Hg-Hg collision, charge 

3 

bohr; C m 

expansion coefficient for n x in Hg-H collision, charge 

3 

bohr; C m 

expansion coefficient for /i x in Hg-Hg collision, charge 

3 

bohr; C m 

expansion coefficient for a in H 0 -H collision, charge 

q x 3 * 

bohr; C m 

expansion coefficient for /j. in Hg-Hg collision, charge 

O 

bohr; C m 

orbital of electron i about proton c, bohr"'^; 
coefficient for quadrupole- induced contribution to x„ com- 

■j 5 

ponent of dipole moment in H 0 -H collision, charge bohr ; 
C m D 

coefficient for quadrupole-induced contribution to x„ com- 

5 

ponent of dipole moment in Hg-Hg collision, charge bohr ; 
C m 5 

expansion coefficient for in H 0 -H collision, charge; C 

ax 3 ^ 

expansion coefficient for in Hg-Hg collision, 

3 

charge; C 
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D m {l = 0, 2, 4) 

T>c &ll 1 l l 2 lJ -2 
d(i)(i= 1,2, 3, 4) 

d z0 (3 = 0, 2, 4) 

h^l^2^2 

dVj(i = 1, 2, 3, 4) 

dvj 

E 

F 

f.(i = 1, 2, 3, 4) 
G a21 

G j 

111 

G q21 


expansion coefficient for n' in H„-H collision, charge; C 

x 3 * 

expansion coefficient for n' in H 0 -H 0 collision, charge; C 

Xg Z Z 

expansion coefficient for n' in Hg-Hg collision, charge; C 

u 

orbital of electron i about proton d, bohr~'* / ^; m" 


coefficient for quadrupole- induced contribution to x Q component of 

£ 4 

derivative of dipole moment in H 9 -H collision, charge bohr ; 

4 z 

C ni 


coefficient for quadrupole- induced contribution to x 0 component of 

O A 

derivative of dipole moment in Hg-Hg collision, charge bohr ; 
Cm 4 

element of volume for electron i in configuration space and spin 
3 3 

space, bohr ; m 

Q *3 

element of volume for electron 1 in configuration space, bohr ; m 
denominator of expression for component of dipole moment 
electric field strength vector, hartree bohr -4 charge -4 ; J m -4 c -4 
constants for fairing of expansion coefficients 

expansion coefficient for overlap contribution to charge bohr, 

C m 

expansion coefficient for charge bohr; C m 

expansion coefficient for quadrupole- induced contribution to M(i)> 
charge bohr; C m 


g 21 

H a21 

H ? 

Ill 

H q21 


coefficient for quadrupole- induced contribution to (1) component of 

5 5 

dipole moment, charge bohr ; C m 
expansion coefficient for overlap contribution to charge; C 

expansion coefficient for charge; C 

expansion coefficient for quadrupole- induced contribution to 
charge; C 

coefficient for quadrupole- induced contribution to (1) component of 

4 4 

derivative of dipole moment, charge bohr ; C m 
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K 

J n fa)(n= 1 , 2 >- 
L (cr,r) 

M(cr) 

M(cr,r) 
m Jaa 0 = 1, 3) 
m jab (i = 1, 3) 

*S 

N l» N 2’ N 3 
Q ’ Q ab’ Q cd 

Q XX’ Q YY’ ^ZZ 

R 

r 

r ab 

r al 

®ab 

s ac 

T(ff) 

T(cr, r) 

U(a) 

U(a, r) 

X,Y,Z 


Xj(j = 1,2,3) 
x n (j = 1, 3) 


Planck constant divided by 2n, J/sec 
. , -5) function of a 

function of a and r 
approximation to M(a, r) for small r 
function of a and r 

integral of square of orbital a(l) and x j a j, bohr; m 
integral of orbitals a(l) and b(l), and x j ab j, bohr; m 
numerator of expression for n for Hg-H, charge bohr; 


C m 


J 


contributions to numerator of expression for n for H 0 -H 9 , 

Xq Cl Cl 

charge bohr; Cm 6 

scalar quadrupole moment of H 2 molecule (subscripts indicate 
which molecule), charge bohr^; C m^ 

0 

elements of H 9 quadrupole moment tensor, charge bohr ; 

? ^ 

C in 

intermolecular distance, bohr; m 

internuclear distance, bohr; m 

distance from proton a to proton b, bohr; m 

distance from proton a to electron 1, bohr; m 

overlap integral for orbitals a(l) and b(l) (in same molecule) 

overlap integral for orbitals a(l) and c(l) (in different mole- 
cules) 

approximation to T(a,r) for small r 
function of o and r 
approximation to U(a,r) for small r 
function of a and t 

Cartesian coordinates with origin at midpoint of line connect- 
ing protons a and b and with Z running along line ab, 
bohr; m 

Cartesian coordinates (see figs. 1 and 5), bohr; m 
Xj coordinate of electron i, bohr; m 
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x 3al 

Xg component of vector from proton a to electron 1, bohr; m 

x 3acl 

Xg component of vector from center of line ac to electron 1, 
bohr; m 

x 3p 

Xg coordinate of proton p, bohr; m 

V 

spherical harmonic 

a(i)(i= 1,2, 3, 4) 
"a 

spin eigenfunction of electron i with component of spin angular 
momentum along axis of quantization equal to K/2 

2 2 -122-1 

polarizability of H atom a, charge bohr hartree ; C m J 

a> a ab’ a cd 

average polarizability of Ho molecule (subscripts indicate which 

-I QOI 

molecule), charge 2 bohr 2 hartree” ; C; nr J” 

/3(i)(i= 1,2, 3,4) 

spin eigenfunction of electron i with component of spin angular 
momentum along axis of quantization equal to -h/2 

A ’ A ab’ A cd 

anisotropy of polarizability of H 0 molecule (subscripts indicate 

2 z 2 -122-1 
which molecule), charge bohr hartree ; C m J 

^a 

orbital exponent of H atom a, bohr” *; m” * 

^ab’^bc’ ^cd 

orbital exponent of H 2 molecule (subscripts indicate which mole- 
cule), bohr"*; m“* 

? 

orbital exponent of two repelling H atoms used in calculating A, 
bohr”*; m”* 

0, 

111 

normalized associated Legendre function 

9, 0 p 0g 

polar angles (see figs. 1 and 5), deg 

^abc 

angle between protons b and c at proton a (first or second quadrant), 
deg 

0 abx.^ = *’ 3 ) 

angle at proton a between proton b and line through proton a parallel 
to the Xj axis (first or second quadrant), deg 

0 abl 

angle between proton b and electron 1 at proton a (first or second 
quadrant), deg 

K > K ab’ K cd 

attractive distortion parameter of H 0 molecule (subscripts indicate 

-11^ 

which molecule), bohr ; m 

X ’ X ab 

repulsive distortion parameter (subscripts indicate repelling orbit- 
als), bohr”*; m”* 


dipole moment vector, charge bohr; C m 
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^a 

M - 1,2,3) 


M 


q x 3 


M x G = 1,2,3) 
X j 


<Pl> ^2 

^cabx^ “ ^ ^ 

* 

^(i = 1, 2, . . . , 5) 
Subscript: 

Hi 

Superscript: 


overlap contribution to /i, charge bohr; C m 
x. component of , charge bohr; C m 

Xg component of quadrupole-induced contribution to ~jl, charge 
bohr; C m 

x^ component of ~ji, charge bohr; C m 

complex components of ~jx, charge bohr; C m 
a with any subscripts 

r times average value of £ for orbitals a(l) and c(l) 

ciC 

t with any subscripts 

dC 

function of r and difference of orbital exponents of orbitals 

clC 

a(l) and c(l) 

azimuthal angles (see figs. 1 and 5), deg 

dihedral angle between plane cab and plane containing a, b, 
and line parallel to the x^ axis, deg 

antisymmetric system wave function, bohr - ; m - 

— _ P. 

determinantal wave function, bohr” ; m - 
configuration i (figs. 2 and 6) 

indicates partial differentiation with respect to internuclear dis- 
tance of left-hand molecule (figs. 1 and 5) 
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